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Abstract — Two-dimensional maps can model inter- 
actions between populations. Despite their simplicity, 
these dynamical systems can show some complex situ- 
ations, as multistability or fractal boundaries between 
basins that lead to remarkable pictures. Some of them 
are shown and explained here for three different 2D 
discrete models. 
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boundary. 

1. Introduction 

Two-dimensional maps can be used to model inter- 
actions between two different species. Such applica- 
tions can be considered in Ecology, Biology or Econ- 
omy [TJ [5] . Generally, real systems consist in a large 
number of interacting species but the understanding of 
the behaviour of such systems in the low dimensional 
case can be of great help as a first step attempt. In this 
work, we consider two-dimensional (2D) models based 
on logistic multiplicative coupling [3 J where complex 
behaviours occur such as multistability phenomena, 
fractal basins of attractors and fractal boundaries be- 
tween basins [4j [5] . These phenomena lead to some re- 
markable graphical representations in the phase space 
plane. In Section 2, we recall three of these considered 
2D models. In Section 3, we emphasize these models, 
which permit to obtain fractal basins. Section 4 is de- 
voted to multistability phenomena and fuzzy or fractal 
boundaries beetwen basins. 

2. The models 

The first considered model is the noninvertible 2D 
map T\ defined by: 



f Xfc+i = A(3a; fc + l)y k {\ - 
\ Vk+i = A(3y fc + l)x fe (l - 
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where A is a real control parameter, x and y are 
real state variables. Previous studies of ([T]) have been 
done in [5]. This model is the symmetrical case of a 
2D model 
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Figure 1: Fractal basins of two order 2 cycles for the 
map T\ ([T]). one basin is in yellow, the other in red. 
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Figure 2: Enlargment of Figure [TJ Zl area. 
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Figure 3: Enlargment of Figure [T] Z2 area. 

proposed for the symbiosis interaction between two 
species [B|. The second model is the noninvertible 2D 
map Ti defined by: 



Xk+i = A(3y fe + l)x k (l - x k ) 
Vk+i = A(3x fc+ i + l)yk(l - Vk) 



(2) 



where A is a real control parameter, x and y are 
real state variables. The map (J2J is also inspired in 
the symbiosis case [6] by including a time asymmetric 
feedback. Previous studies of © have been presented 

in0. 



.T fc+ i = A(-3y fc + 4)x fc (l - a; fc ) 
2/fc+i = A(-3a;A; + i)y k (l - y k ) 



(3) 



where A is real, x and y are real state variables. 
The map ((3]) corresponds to a competitive interaction 
between two species [8]. In each model, A measures 
the strength of the coupling. 
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Figure 5: Basin of a fixed point for the map T2 (J2|). 
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Figure 4: Enlargment of Figure Z21 area. 
The third model is the noninvertible 2D map T3 



Figure 6: Enlargment of Figure [SJ Fractal structure is 
observed. 



3. Fractal basins 

Figures [T]|4] show the basin of two coexisting attrac- 
tors and successive enlargements in the case of the 
map Ti, each attractor being an order 2 cycle: one 
basin is in yellow, the other one in red. Basins are 
symmetrical, non-connected and fractal. Successive 
enlargements of the 2-dimensional phase space on the 
diagonal (Zl area) and on Z2 area show auto-similarity 
properties. The shape of each enlarged piece of basin 
is similar to the precedent piece with an alternation 
between the successive locations of symmetrical yel- 
low and red basins. Such a shape can be explained by 
using the critical manifolds of (TTJ) [3J [3] . 

Figures [5][6] show the fractal basin of a fixed point 
for the map Ti- This basin is non-connected. Figure|5] 
shows clearly the auto-similarity property. The ap- 
pearance of such a fractal basin can also be explained 
by using the critical manifolds [4J [Jj . 

Figures [7][8] show a fractal basin for the map T3. 
The attractor is chaotic and it is not represented on 
the Figures [8]. As in the case of the map ([T]), the 
successive auto-similar pieces of the basin are located 
on the diagonal. The fractalization of the basins can 
be understood by means of the critical manifolds |4j [8] , 
as in the case of the maps (p} and ((2|). 
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Figure 7: Fractal basin for the map T3 ([3|). 



4. Multistability and Fractal boundaries 

Multistability among several attractors is very fre- 
quent in these models. As it can be seen in the Fig- 
ures l9"][T2l the basins can take complex and fractal 
forms. Figure [5] shows the shape of the basins of three 
different attractors. The boundary among those basins 



is fractal, due to the accumulation and the entangle- 
ment close to the boundary of the square [0, l] 2 . Fig- 
ure [TU] represents the basin of two different attractors, 
wich are order 3 invariant closed curve (ICC), for the 
map T-2. Then we obtain six different basins in T% 
for each piece of the two co-existing order 3 ICC in 
Ti- Some of the basins can be riddled, that is, suc- 
cessive zooms of a basin zone cannot differentiate the 
boundaries between the basins of the different attrac- 
tors. Such basins are shown on Figures ITTl[T2l There 
are two chaotic attractors, which are order 52 chaotic 
rings. They are called weakly chaotic because of the 
greater Lyapunov exponent, which is slightly positive 
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Figure 8: Enlargment of Figure [71 
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Figure 9: Basins of an order 29 cycle and two order 3 
invariant closed curves (ICC) for the map T 2 . 
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Figure 10: Basin for T| of 2 order 3 ICC. 
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Figure 11: Riddled basin of two order 52 weakly 
chaotic rings (WCR) for the map T3. 
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Figure 12: Enlargment of Figure [TTJ BZ area. 
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